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The statistical properties (fluctuations) of heartbeat intervals (RR intervals) in 
ECG are studied and compared with the predictions of Random Matrix Theory 
(RMT). It is found that heartbeat intervals only locally exhibit the fluctuation pat- 
terns (universality) predicted by the RMT. This finding shows that heartbeat dy- 
namics is of the mixed type where regular and irregular (chaotic) regimes coexist and 
the Berry-Robnik theory can be applied. It is also observed that the distribution of 
heartbeat intervals is well described by the one-parameter Brody distribution. The 
parameter (3 of the Brody distribution is seen to be connected with the dynamical 
state of the heart. 



1 Introduction 

The time series of heartbeat intervals (RR intervals), used in various analy- 
ses, are usually Holter type data or data from steady state ambulatory mea- 
surements. The great wealth of data about the dynamics of the heart that 
is contained in such ECG records is usually reduced to characterize only the 
mean heart rate and the presence and frequency of some abnormal electrocar- 
diographic complexes. The analysis of short-range and long-range fluctuations 
and its universality are largely ignored. ^ The normal human heartbeats display 
considerable variability and nonstationarity over time. The result is that fluctu- 
ations of heartbeat intervals around the mean value are irregular and complex, 
even under resting conditions, Fig. 1. A number of new methods have been 
proposed to quantify these complex fluctuations, mainly based on the nonlinear 
and fractal analysis. 2 )> 3 )' 4 )' 5 ) 

In this paper we examine the possibility that the distribution of heartbeat 
intervals exhibits universality of the type predicted by the RMT. 6 ) 
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2 Statistical analysis of heartbeat intervals 



The heartbeat time series of an human represents an ordered sequence of R- 
waves (QRS-complexes) occuring at times {ii, £2, • • • , U, • • •}. In this form the 
heartbeat dynamics can be regarded as a fractal stochastic point process on a 
time-line. The point process is completely described by the set {t,}, or equiva- 
lently by the set of RR-intervals defined as {RR(i) = t i+1 — U, % — 1, 2, • • •}. 

To see the universality in beat-to-beat RR interval fluctuations we need to 
magnify the heartbeat spectrum {tj} so that mean spacing is unity. Let us first 
determine the number of heartbeats up to time t. It is given by the staircase 
function N(t): 

N{t) = "£e(t-u) (i) 

i 

where 9(t) is the step function. 

In order to be able to compare statistical properties of heartbeat spectra, 
in particular the fluctuations around the mean value, with the predictions of 
the RMT it is necesary to normalize (unfolde) the time spectrum {£;} to {71} 
in such a way that the mean spacing between neighbouring r« is unity. This is 
achieved by the mapping 

Ti = N(ti), Si = T i+1 -n (2) 

where N(t) denotes the smoothed staircase which account for the mean trend 
in N(t) and Si is the normalized heartbeat spacing. The "local" mean spacing 
is defined as A(£) = [dN (t) / dt]" 1 so that Si becomes 

t i+ i - U 

1 A[(* i+ i +U)/2]' [6) 

During the real time measurements, we always deal with the finite number of 
heartbeats, say n. In that case, we define the "global" mean spacing as 

1 n ~ l f — + 

A = _i_g M(0 = ^l (4) 

so that St = RR(i)/A. 

There are two important statistical measures that we shall use to charac- 
terize the heartbeat spectra {rj}. One is the nearest-neighbor heartbeat spac- 
ing distribution P(S)dS which gives the probability to find in the interval 
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(S, S + dS). P{S) is a short range statistics, it measures the short range corre- 
lations between heartbeats. Its behaviour for small spacing, S — > 0, 

P(S) « const x ^ (5) 

determines the type of short range correlations between heartbeat intervals. 
Two cases can be distinguished: a) when (3 = meaning that the short range 
fluctuations of {Si} are just those of uncorrelated random numbers, i.e. Pois- 
sonian with a well-known phenomenom of attraction (clustering), and b) when 
(3 7^ meaning that the long range fractal fluctuations are present in the se- 
quence of Si that is characterized by the phenomenom of repulsion. The uni- 
versal statistics (no free parameters) predicted by the RMT 6 ) are given by 

Ppoisson(S) = exp(-S), (3 = 0, 
Pgoe(S) = ^Sexp(~S 2 ), j3 = l, (6) 

Pgue(S) = ^S 2 exp(--S 2 ), (3 = 2, 

7T 7T 

The GOE (Gaussian orthogonal ensemble) distribution, known as a Wigner 
distribution, should be valid for irregular (chaotic) systems with an antiunitary 
symmetry, whereas GUE (Gaussian Unitarian Ensemble) should apply if the 
system has no such symmetry. 

To avoid the smoothing procedure on a measured data, it is easier and less 
ambiguous to study the inegrated or the cumulative RR intervals 

I(S)= [ S p(S)dS (7) 
Jo 

which yields a useful statistics even for small sample of RR intervals. This gives 

Ipoi S son(S) = 1 - eXp(-S), Iwigner(S) = 1 - 6Xp(--S 2 ), (8) 

and a nonanalytical expression for Igue(S). 

Another important statistics, not considered here, is the spectral rigidity 
A3 of Dyson and Metha 7 - ) which characterizes the long-range correlations of 
heartbeats. It can be related to the Hurst exponent 8 - ) ' 9 ^ ) associated with the 
heartbeat time series. 



3 




Figure 1: Time series of RR intervals as a function of beat number from a 
person during an ergometric measurement. 

3 Results and discussion 

Time series of RR intervals were measured at the Institute for Cardiovascular 
Disease and Rehabilitation in Zagreb. Each masurement of ECG was over a 
time duration of about 15 min. (~ 2000 heartbeats), in a controlled ergometric 
regime. These type of measurements are used as routine in everyday clinical 
diagnostic practice, because some heart deseases, such as stable angine pectoris 
(SAP), usually become transparent under physical activities. 

The ECG data were digitized by the WaveBook 512 (Iotech. Cal. USA), 
and transferred to a computer. The RR interval series was passed through a 
filter that eliminates noise and artefacts. All R-wave peaks were first edited 
automatically, after which a careful manual editing was performed by visual 
inspection of the each RR interval. After this, all questionable portions were 
excluded manually, and only segments with > 90% sinus beats were included in 
the final analysis. 

Each measurement consisted of stationary state part (pretrigger Pt), a few 
stages of running (P1-P4) and a period of relaxation (Re), Fig. 1. Here we 
report on the universality analysis of RR interval fluctuations in the pretrigger 
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and PI period only. 10 ) Our patients were divided in two groups: one with 
the evidence of ishemic ST-segment depression (SAP subjects), and the control 
group of healthy subjects. Selection of subjects was performed by a cardiologist 
according to the generally accepted medical knowledge. 

The integrated RR interval distributions I(S) are shown in Fig.2 together 
with Iw{S). Clearly the universal, parameter free, distributions Pp isson,Pw 
and Pque of RMT are excluded. A possible reason for this is that heart is a 
very complex system of nonlinearly coupled biological oscillators whose motion 
is partly regular and partly irregular (chaotic) or mixed . Another possibility 
is that universality is only local. 10 ) 

The first theoretical step towards the understanding of the regular-irregular 
coexistence in a dynamical system was the work by Berry and Robnik ll > in 
which they offered a semiclassical model for P(S) known as the Berry- Robnik 
formulae. It is to be applied at large intervals S > 1. At small and intermediate 
intervals S < 1 the Brody distribution 12 ) 

P B (S) = aS^exp(-bS^ 1 ), a = (J3 + l)b, b = {r(|±|)^ +1 (9) 

gives better fitting results. The Brody distribution becomes Poisson at (3 — 
and Wigner for {3 = 1. 

Figure 3 shows that the Brody distribution fits very well the ECG measured 
distribution of RR intervals. Each subject has its own value of (3 which might 
be an indicator of a dynamical state of the heart. Indeed, Fig. 4 clearly shows 
that (3 could be a useful parameter for separating SAP subjects from healthy 
ones. 

4 Conclusion 

In this paper we have investigated the possibility that fluctuations of RR inter- 
vals exhibit universal statistical laws. It is found, Fig. 2, that this is not so. 
Heart is a complex nonlinear system with many coexisting regular and irreg- 
ular motions which manifest themselves through the nonuniversal fluctuations 
of RR intervals. It is also found, Fig. 3 and 4, that the one-parameter Brody 
distribution could be used successfuly to describe the fluctuation pattern of 
RR intervals. The parameter (3 of the Brody distribution is seen, Fig. 4, to be 
connected with the dynamical state of the heart. 

Further studies in larger populations are needed to confirm these interesting 
results. 
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Figure 2: Integrated RR interval distributions I(S), in PI period, in comparison 
with Wigner distribution. 



6 






Figure 3: RR interval distributions, in Pt period, from three different per- 
sons in comparison with the one parameter Brody distribution. [S = (RR — 
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Figure 4: Brody parameter (3 as a function of a healthy state of the heart. 
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